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Abstract. We report the experimental observation of the disruption of the superfluid 
atomic current flowing through an array of weakly linked Bose-Einstein condensates. 
The condensates are trapped in an optical lattice superimposed on a harmonic magnetic 
potential. The dynamical response of the system to a change of the magnetic potential 
minimum along the optical lattice axis goes from a coherent oscillation (superfluid 
regime) to a localization of the condensates in the harmonic trap ( "classical" insulator 
regime) . The localization occurs when the initial displacement is larger than a critical 
value or, equivalently, when the velocity of the wavepacket's center of mass is larger 
than a critical velocity dependent on the tunnelling rate between adjacent sites. 
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Atomic Bose-Einstein condensates have been either loaded or produced in periodic 
potentials opening up the possibility to investigate new phenomena tuning the degree of 
coherence in the system. Experiments have explored regimes ranging from the coherent 
matter wave emission from a condensate loaded on a vertical standing wave pQ, to the 
observation of number squeezed states [2], to the demonstration of a one-dimensional 
Josephson junctions array with a linear chain of condensates produced in an optical 
lattice [3], to the recent observation of a quantum phase transition in a condensate loaded 
in a 3D optical lattice [3]. The dynamical behavior of coherent matter waves in periodic 
potentials has also been the subject of extensive theoretical work. Different phenomena 
have been predicted including the formation of bright solitons and the appearance of 
various forms of instabilities where the usual long range phase coherence properties are 
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In this paper we report on the experimental observation of the disruption of 
superfluidity in a linear array of 87 Rb condensates trapped in the combined potential of 
an optical standing wave superimposed on a harmonic magnetic trap. 

We observe a transition from a regime in which the wavepacket coherently oscillates 
in the array to another one in which the condensates stop in the optical potential 
sites and lose their relative phase coherence. We also observe the appearance of an 
intermediate regime where no oscillations are visible while the relative phase coherence 
is not completely lost. In a previous experiment we have investigated the regime of small 
periodic potential height, showing that in this regime and for large atom numbers, the 
mechanism for superfluidity breakdown was in agreement with predictions based on 
the onset of the Landau instability [T5] . In the present experiment we are exploring 
the dynamics of the system in the case of a large periodic potential height, where a 
modulational instability is expected to play a fundamental role [HJ Ej. Furthermore, 
one can study the overall coherence of the system observing the phase evolution in the 
interferogram of the expanded array of condensates. 

The experimental setup has been already described in a dilute vapor of 87 Rb 
atoms confined in a Ioffe-type magnetic trap in the (F = 1, rriF = —1) state is Bose 
condensed by the rf-evaporative cooling. Slightly above the condensation threshold, 
we superimpose to the harmonic magnetic trap a ID optical standing wave. We 
then continue the evaporation process through the phase-transition temperature. The 
magnetic trap is cigar-shaped, with frequencies uo z = 2tt x 9 Hz and u r = 2tt x 90 Hz along 
the axial and radial directions, respectively. The optical standing wave is superimposed 
to the magnetic potential, along the z-axis, and is created with a retroreflected 
collimated laser beam, detuned ~ 3 nm to the blue of the rubidium Dl transition 
at A = 795 nm. The resulting periodic potential is Vl = Vocos 2 (2ttz/\): the valleys 
of the potential are separated by A/2 and the interwell energy barrier Vq — sEr can 
be controlled by varying the intensity of the laser beam. Er = h 2 /2m\ 2 is the kinetic 
energy of an atom of mass m recoiling after absorbing one lattice photon. Typically, 
in our experiments, s ranges from 3 to 14 Er, for such values, the condensate chemical 
potential //, that ranges from 1.8 to 3.7 E R , is smaller than V and the system realizes 
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an array of weakly coupled condensates driven by an external harmonic field. 

We drive the system out of equilibrium by a sudden displacement (t dis <C 2n/oj z ) 
of the magnetic potential along the lattice axis. The ensuing dynamics is revealed by 
turning off, after different evolution times, both the magnetic and optical traps and 
imaging the atomic density distribution after an expansion time t exp = 27.8 ms. During 
the expansion the different condensates, originally located at the nodes of the laser 
standing wave, overlap and interfere allowing to monitor their relative phases [TT)] . 

For small displacements we observe a superfluid regime that has been extensively 
studied both at zero temperature [3] and in the presence of a significant thermal fraction 
[T7j . The condensates released from the traps show an interference pattern, consisting of 
three peaks, separated by a distance d = 2ht exp /m\, oscillating in phase, thus revealing 
the long-range coherence of the condensates across the entire optical lattice. 

When we increase the displacement of the magnetic trap center for a fixed height 
of the optical barriers we observe a different behaviour: the center of mass of the atomic 
sample no longer oscillates but slowly moves towards the center of the magnetic potential 
where it subsequently stops. This is shown in Fig.^where we report the evolution of the 
center-of-mass position of the atomic cloud inside the combined trap for a barrier height 
of 5 E# and displacements of 30 /im and 120 /im. We reconstruct the center of mass 
motion in the trap from the images of the expanded cloud assuming ballistic expansion. 
We take data points at intervals of time small enough so as to assume the condensate 
moving in the combined trap with constant velocity during each time interval. 
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Figure 1. Position of the center-of-mass of the atomic cloud in the trap for an optical 
barrier height of 5 E/?. Filled circles: displacement of 30 /im; open circles: displacement 
of 120 /im. The continuous line represents a fit to the "small displacement" data with 
a sine function. 
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For the larger displacement we also observe a vanishing of the visibility of the 
interference peaks demonstrating the loss of long-range coherence across the array. 
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Figure 2. Axial width (R z ) of the interferogram central peak normalized to the peaks 
separation d as a function of time spent in the combined trap after the magnetic trap 
displacement for an optical barrier height of 5 E#. Filled circles: displacement of 
30 /im; open circles: displacement of 120 /im. In the right part of the figure we show 
the observed interferogram at t=120 ms for the two displacements. 

In Fig. |21 we report the evolution of the axial width R z of the interferogram central 
peak normalized to the peaks separation d for the same conditions as in Fig. ^ While 
in the small amplitude coherent oscillation the axial width remains constant, in the 
larger displacement case it spreads out reaching the peak separation in ~120 ms when 
the interferogram is completely washed out. We remark that the vanishing of the 
interferogram is a dynamical effect and is not due to scattering of lattice photons. 
Indeed in the absence of displacement or even for small displacements we do not observe 
a deterioration of the interference contrast on the time scales of the experiment as shown 
in the right part of Fig. El 

We have studied the crossover between this two dynamical regimes repeating the 
experiment varying the displacement for different optical lattice heights Vq. In Fig. El 
we summarize the experimental results: the closed circles correspond to values of 
displacement and lattice depth for which we observe a coherent oscillation. The open 
circles identify conditions where we have the destruction of the oscillation and a slow, 
overdamped, motion toward the center of the magnetic trap. In the latter case, only for 
very large displacements we observe a complete vanishing of the interferogram. 

We compare our results with the prediction for the onset of the dynamical instability 
recently proposed in [Hj using the one dimensional Discrete Non-Linear Schrodinger 
equation (DNLS) model [3]. This treatment gives an analytical expression for the 
"critical" velocity v cr to enter the dynamical instability regime. The ID system is 
predicted to undergo a sudden transition from a regime with coherent oscillations 



5 



(superfluid regime) to one with pinning (insulator) when the phase difference between 
adjacent condensates A0(i) reaches it/2. In this case the relative phases start to run 
independently with different velocities, the phase coherence through the array is lost 
and the system behaves as an insulator. The system reaches this regime for a critical 
velocity 

«c=- r - (1) 

where K is proportional to the tunnelling rate between adjacent sites of the optical 
potential. In a harmonic potential this corresponds to a critical displacement: 



A, cr . = -^- (2) 

where Q = ^muj^(^) 2 describes the magnetic potential energy of the condensates. Since 
the tunnelling rate K depends on Vb, the critical displacement should depend on the 
height of the interwell potential: the greater Vq is, the lower the tunnelling rate and the 
critical displacement are. 

In Fig. El we report the DNLS prediction for the critical displacement (continuous 
line JBj)- Effectively, the theoretical curve well divides the two dynamical regions 
observed in the experiment. 




Figure 3. Comparison of experimental results with Eq. ©. Filled circles represent 
coherent oscillations (with coherent interferogram featuring the three peaks), empty 
circles denote pinned motions; continuous line results of Eq. @. 



From the experimental data (open circles in Fig. |3J) we can also extract the 
maximum velocity reached by the system before entering the instability regime. In 
Fig. HJ we show the experimental maximum velocities measured as a function of the 
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optical potential height, in qualitative agreement with the theoretical prediction of 
Eq. (HJ) (continuous line). 
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Figure 4. Critical velocity as a function of the optical lattice depth. The experimental 
results (filled circles) are compared with with Eq. Q (continuous line) . 

We believe that the ID model gives a reliable estimate for the onset of the loss 
of superfluidity, however it does not quantitatively describe the degradation of the 
interferogram in the present experimental setup. Actually, in the experiment, for values 
of displacement just above the predicted Az cr we could not observe a complete disruption 
of the interferogram but rather the appearance of more complex structures within the 
three visible peaks. We believe that this could be due to the excitation of radial modes 
that are obviously not included in the simple ID model. A possible way to include such 
modes is to resort to a full numerical solution of the 3D Gross-Pitaevskii equation |19j . 

In conclusion, we have experimentally observed the transition between different 
regimes in an array of Josephson junctions realized with BECs trapped in a one 
dimensional optical lattice. The transition occurs between a superfluid and an insulator 
regime and is accompanied by a loss of coherence through the array even though each 
condensate in the array is still described by a coherent state. The experimental findings 
have been compared with the prediction of a one dimensional theoretical model [B] based 
on a discrete nonlinear Schrodinger equation where the transition is due to a dynamical 
instability taking place when the eigenfrequencies of the excitation spectrum become 
imaginary. This model is found to qualitatively describe the instability onset which 
occurs at a critical value of the displacement or, equivalently, when the velocity of the 
wavepacket's center of mass is larger than a critical velocity dependent on the tunnelling 
rate through the optical barriers. We have also observed that close to the transition the 
ID model fails to describe the experimental results where the system behaves like an 
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insulator but coherence is still present through the array. We attribute this behavior to 
the excitation of radial modes that would require a full 3D simulation of the GPE. 

A quantitative comparison between the critical velocity for the dynamical instability 
and the speed of sound in the presence of the optical lattice in the regime of high optical 
lattice to our knowledge is still lacking. We think that it will be very helpful to enlighten 
the role of this two instability regimes already discussed for the case of small lattice 
height in 0. 

We are indebted to F. Minardi and S. Burger for their help in the initial stage of 
the experiment and for several discussions. We acknowledge A. Smerzi, A. Trombettoni, 
A. R. Bishop, P. G. Kevrekidis, M. Modugno, A. Tuchman and F. T. Arecchi for many 
fruitful discussions. Experimental evidences of a similar instability has also been found 
in j2B|. 
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